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Novel analysis of finite dimensional Hilbert space is outlined. The approach bypasses general, 
inherent, difficulties present in handling angular variables in finite dimensional problems: The finite 
dimensional, d, Hilbert space operators are underpinned with finite geometry which provide intuitive 
perspective to the physical operators. The analysis emphasizes a central role for projectors of mutual 
unbiased bases (MUB) states, extending thereby their use in finite dimensional quantum mechanics 
studies. Interrelation among the Hilbert space operators revealed via their (finite) dual affme plane 
geometry (DAPG) underpinning are displayed and utilized in formulating the finite dimensional 
ubiquitous Radon transformation and its inverse illustrating phase space-like physics encoded in 
lines and points of the geometry. The finite geometry required for our study is outlined. 

o 

I. INTRODUCTION 

(D, 

[j^ . The mathematics of Radon transform was introduced for studies of Astronomical data [l], |2g, [35| . Today it finds an 
extensive application in CAT (computer assisted tomography) scan ,state reconstruction [!, I8l-fl0j as well 

as within the fundamental phase space formulation of quantum mechanics |3l. l3l| . The methodology is geometrically 
based: A phase space image of an arbitrary Hilbert space operator A is given by 33] , 
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^— > 



A -> W A (q,p) = J dye^iq - y/2\A\q + y/2). (1) 

The Wigner function (32|, W p (q,p), being the image of the density operator, p, although not non- negative, does 
enjoy many of the attributes of a phase space distribution. In particular its marginal, 



p{x',9) = dqdpS(Cq + Sp-x')W p (q,p), C = cos6,S = sin9, (2) 



> 

^2 ' is its Radon transform [3, 23]: 

n[W p ](x,6)=p(x,6); K- 1 lp}(q,p) = W p (q,p). (3) 

The explicit expression for W p (q,p) in terms of principal value integral of p(x,9 is given in pit [8HTfl[35j. We offer 
04 ■ in the last section an interpretation for this singular relation gained via the finite dimensional approach. 



The tomographic luster of the Radon transform suggests viewing it via mutual unbiased bases (MUB). Informa- 
tionally complete set of MUB, i.e. such that allows complete accounting of an arbitrary state, are {|a;',0)}; — oo < 
x' < oo, ; < 9 < it. These are the eigenstates of the dynamical variables Xg — Cx + Sp, [13, [35[. Now 

trpmathbbP(x', 9) = J dqdpS(Cq + Sp - x')W p {q,p) = p(x\ 9), (4) 

illustrating a quasi distributional attribute of W p (q,p), as the phase space mappings of F(x',9) = \x\ 9) (9, x'\ is 
S(x' — Cq — Sp). Noteworthy is the pivotal role of the delta function: it identifies the Radon transformation, defines 
the marginal distribution and constitute the phase space mapping of an MUB projector. 

The present intuitive approach involves the replacement of phase space with finite geometry's points and lines. 
The role of the Dirac delta function is played by what we signify with A a _j where a and j refer to finite geometry 
point and line respectively. In this way we eschew in the finite dimensional intrinsically sticky issue of angular 
variables. On the physical side the salient feature is the mapping of finite dimensional Hilbert space operators 
onto c number functions of the phase space like geometrical points and lines and, in particular, the definition of 
quasi distribution and Radon transformation over this phase space like points and lines. These ideas are clarified below. 
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II. FINITE GEOMETRY AND HILBERT SPACE OPERATORS 

We now briefly review the essential features of finite geometry required for our study 0, mm ©si- 

A finite plane geometry is a system possessing a finite number of points and lines. There are two kinds of finite 
plane geometry: affine and projective. We shall confine ourselves to affine plane geometry (APG) which is defined 
as follows. An APG is a non empty set whose elements are called points. These are grouped in subsets called lines 
subject to: 

1. Given any two distinct points there is exactly one line containing both. 

2. Given a line L and a point S not in L (S 3 L), there exists exactly one line L' containing S such that Lf]L' = 0. 
This is the parallel postulate. 

3. There are 3 points that are not collinear. 

It can be shown [28l - l30l | that for d — p m (a power of prime) APG can be constructed (our study here is for d=p) and 
the following properties are, necessarily, built in: 

a. The number of points is d 2 ; S a , a — l,2,...<i 2 and the number of lines is d(d+l); Lj, j — 1, 2....d(d + 1). 

b. A pair of lines may have at most one point in common. 

c. Each line is made of d points and each point is common to d+1 lines: Lj = \J a S^, S a — C(j=\ ■ 

d. If a line Lj is parallel to the distinct lines Lk and Li then Lk \\ Li. The d 2 points are grouped in sets of d parallel 
lines. There are d+1 such groupings. 

e. Each line in a set of parallel lines intersect each line of any other set. 



The existence of APG implies |21l l28H3G |the existence of its dual geometry DAPG wherein the points and lines are 
interchanged. Since we shall study extensively this, DAPG, we list the corresponding properties for it. We shall refer 
to these by DAPG(y): 

a. The number of lines is d 2 , Lj, j = 1, 2....d 2 . The number of points is d(d+l), S a , a — 1,2, ...d(d + 1). 

b. A pair of points on a line determine a line uniquely. Two (distinct) lines share one and only one point. 

c. Each point is common to d lines. Each line contain d+1 points. 

d. The d(d+l) points may be grouped in sets, R a , of d points each no two of a set share a line. Such a set is 
designated by a' £ {a U M a }, a' — 1, 2, ...d. (M a contain all the points not connected to a - they are not connected 
among themselves.) i.e. such a set contain d disjointed (among themselves) points. There are d+1 such sets: 

d(d+l) d 

\J S a =\jR a ; R a = |J S a ,; R a f)R a ,=0,a^a'. (5) 

a=l a=l a'ea\JM a 

e. Each point of a set of disjoint points is connected to every other point not in its set. 
DAPG(c) allows the identification, which we adopt, 

Sa=~X>. =* E S a , = \Y,Lj. (6) 
We list now some direct DAPG implied interrelation subjected to Eq. ([6]): DAPG(d) implies 

d+i d d(d+i) d 2 

E E s«= E ** = ^E^ =* 

a a'eaUM a 

d(d+l 



E ^ = ^E^ = ^iE^ (?) 



a'EaUM, 



III. FINITE DIMENSIONAL MUTUAL UNBIASED BASES, MUB, BRIEF REVIEW 



In a finite, d-dimensional, Hilbert s pac e two complete, orthonormal vectorial bases, Bi, B2, are said to be MUB if 
and only if (i?x ^ S 2 ) El MM M M M, HE 113 : 
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V|u), \v) e B u B 2 resp., \(u\v}\ = l/Vd. (8) 

The physical meaning of this is that knowledge that a system is in a particular state in one basis implies complete 
ignorance of its state in the other basis. 

Ivanovic [1 31 ] proved that there are at most d+1 MUB, pairwise, in a d-dimensional Hilbert space and gave an explicit 
formulae for the d+1 bases in the case of d=p (prime number). Wootters and Fields [f| constructed such d+1 bases 
for d = P m with m an integer. Variety of methods for construction of the d+1 bases for d = p m are now available 
EH El- Our present study is confined to d = p ^ 2. 

We now give explicitly the MUB states in conjunction with the algebraically complete operators pi [l2| set: Z,X. 
Thus we label the d distinct states spanning the Hilbert space, termed the computational basis, by |n), n — 
0,l,..d- l;|n + d) = \n) 

Z\n) = u/»; X\n) = |n + 1), u = e l2lT / d . (9) 
The d states in each of the d+1 MUB bases [HI, [3]are the states of the computational basis (CB) and 

1 d-l 

|m;t) = T y^"("- 1 )-™|n); 6, m = 0, 1, ..d - 1. (10) 
Vd 5 

Here b labels the other d bases and m labels the states within a basis. We have [l4[ 

XZ b \m;b) =io m \m-b). (11) 

For later reference we shall refer to the computational basis (CB) by b=-l. Thus the above gives d+1 bases, b=- 
1,0,1, ...d-l. The total number of states d(d+l) are grouped in d+1 sets each of d states. We have of course, 

(m;b\m';b) = S m , m ,; |(m; b\m'; b')\ = b^b'. (12) 

We remark at this junction that the eigen values of the CB might be considered finite dimensional modulated position 
values ("q") and the eigenvalues of shifting operator, X, modulated momentum ("p"). 
This completes our discussion of MUB. 



IV. DAPG UNDERPINNING OF D-DIMENSIONAL HILBERT SPACE 



Now the underpinning of Hilbert space o per ators with DAPG will be undertaken. We consider d=p, a prime, + 1 2. 
For d=p we may construct d+1 MUB [1, E3 • Points will be associated with MUB state projectors. To this end 

we recall that we designate the MUB states by |m, b). with b — 0, 1, 2...d — 1 labels the eigenfunction of, resp. XZ b . 
m labels the state within a basis. We designate the computational basis, CB, by b=-l.( Note that the first column 
label of -1 is for convenience and does not designate negative value of a number.) The projection operator defined by, 

A a = \m,b)(b,m\; a = {b,m}; b = -1,0, 1, 2...d- 1; m = 0, 1, 2, ..d - 1. (13) 

The point label, a = (m, b) is now associated with the projection operator, A a . We now consider a realization, 
possible for d=p, a prime , of a d dimensional DAPG, as points marked on a rectangular whose horizontal width 
(x-axis) is made of d+1 columns of points. Each column is labelled by b, and its vertical height (y axis) is made of d 
points each marked with m. The total number of points is d(d+l) - there are d points in each of the d+1 columns. 
We associate the d points m — 0,1,2, ...d — 1 in each set, labelled by b, (a ~ (m, b)) to the disjointed points of 
DAPG(d), R a viz. for fixed ba'EaU M a form a column. The columns are arranged according to their basis label, 
b. The first, i.e. left most, being b=-l, a_i = (m, —1); m = 0, 1, ...d — 1. Next we place the columns with increasing 
values of b, the basis label, as we move to the right. Thus the right most column is for b=d-l. Lines are now made of 
d+1 points, each of different b (necessarily , DAPG(b)). A point S a underpins a Hilbert space state projector, A a . 
i.e. A 2 a — A a , and trA a — 1. We designate the line operator underpinned with Lj, by Pj. Thus the above relations 
now hold with S a O A a ; Lj O Pj . 
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e.g. for d=3 the underpinning's schematics is, 
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Now DAPG(c) (and Eq . (fT3")l . (|7|) ) implies that A a ; a = 0, 1, 2...d — 1; a G a' U M a > forms an orthonormal basis for 
the d-dimensional Hilbert space: 

d d d(d+i) 

J2\™,b){b,m\= A a = {d + l)I. (14) 

m a'eaUMo, a 

Eq.© implies, 



A « = ^E>- ( 15 ) 

Evaluating 



<2 



E^^EE^-^t E Pf + (d+i)Pi] =i+p 3 ^ p, = (i6) 



d ^— ^ ~— ' d 



Eq.jT2]) implies, 



Hence, using Eq.(j6l).(fT2|). 



Trivially 



1; a = a' 

trA a A a , = { 0; a ^ a'; a € a' U M a >. (17) 
j ; a ^ a' ; a a' U M a ' . 



[LaV" trA a A a' - A a = 0; a 3 J. 



d+1 

trPj = Y trA <* -1 = 1- (19) 



d+l 



trP J = E trAa -1 = 1- =>■ fr^i-Pj' = E *r-Pj^a' - 1 = \ „ J , { ( 2 °) 
i.e. 

(21) 

An alternative view of the Lambda function is gained via 



' ' ~ d ^ ' ' " \V /•,•/', = 0; ^ a. 



Note that the case of j : ^ a implies j G M a 
These are summarized by 



(22) 



A Q , . trA aPj = J ; a e J = ^ j ^ Q (23) 
0; «G\?, I 0; j ^ a. 
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V. GEOMETRIC UNDERPINNING OF MUB QUANTUM OPERATORS: THE LINE OPERATOR 

We now consider a particular realization of DAPG of dimensionality d = p, =/= 2 which is the basis of our present 
study. Thus a — m(b) designate a point by its row, m, and its column, b; when b is allowed to vary - it designate 
the point's row position in every column. We now assert that the d+1 points, m,j(b), b = 0, 1, 2, ...d— 1, and rrij(— 1), 
form the line j which contain the two (specific) points m(-l) and m(0). The line is given by (we forfeit the subscript 
j - it is implicit), 

m(b) = -(c- l) + m(0), mod[d\ b^-1, 
m(-l) = c/2. (24) 

The rationale for this particular form will be clarified below. Thus a line j is parameterized fully by j = 
(m(— 1), m(0)). (Note: since our line labelling is based on b values -1 and a more economic label for j is j — (m_i,mo) 
i.e. the m values for b=-l and 0. We shall use either when no confusion should arise.) We now prove that the set 
j = 1, 2, 3...d 2 lines covered by Eq.flU]) form a DAPG. 

1. Since each of the parameters, m(-l) and m(0), can have d values the number of lines is d 2 ; the number of points 
in a line is evidently d+1: a point for each b. DAPG(a). 

2. The linearity of the equation precludes having two points with a common value of b on the same line, DAPG(d). 
Now consider two points on a given line, m(b\), mfa); b\ ^ bi. We have from Eq. (|24")) . (b ^ — 1, b\ ^ 62) 

m(bi) — -^-(c — 1) + m(0), mod[d] 
62 

m{b 2 ) = y(c-l) + m(0), mod[d}. (25) 
These two equation determine uniquely (for d—p, prime) m(-l) and m(0). DAPG(b). 

For fixed point, m(b), c m(0) i.e the number of free parameters is d (the number of points on a fixed column). 
Thus each point is common to d lines. That the line contain d+1 is obvious. DAPG(c). 

3. As is argued in 2 above no line contain two points in the same column (i.e. with equal b). Thus the d points, 
a, in a column form a set R a = Ua'eauM ^ a '> wl ^h trivially R a f]R a i = 0, a / a', and (Ja=i~^ &a = Ua=i 
DAPG(d). 

4. Consider two arbitrary points not in the same set, R a defined above: m(bi), 771(62) (bi ^ 62)- The argument of 
2 above states that, for d=p 1 there is a unique solution for the two parameters that specify the line containing these 
points. DAPG(e). 

We illustrate the above for d=3, where we explicitly specify the points contained in the line j = (m(— 1) = 
(l,-l),m(Q) = (2,0)) 

/ m\b -1 1 2 \ 

(0,2) 

1 (1,-1) ■ (1,1) ■ 

V 2 • (2,0) • • / 

For example the point m(l) is gotten from 

m(l) = ~(2 - 1) + 2 = 1 mod[3] -> m(l) = (1, 1). 
Similar calculation gives the other point: m(2)=(0,2). i.e. the line j=(l,2) contains the points (1,-1), (2,0), (1,1) and 
The geometrical line, Lj, j = (1, 2) given above upon being transcribed to its operator formula is via Eq. (|16l) . 



Pj=(h2) = A (h-1) + A (2,0) + 4(1,1) + 4(0,2) - I- (26) 

Evaluating the point operators, A a , 

fo 0\ 

4(1,-1) =010 ,A ( 2, 0) =5[w 1 u 2 \ ,A (1A) = ^{oj 2 1 1 j ,A (0i 2) = ± I 1 1 w ] , (27) 
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and evaluating the sum, Eq.([26|). gives 



w> 

fj:(m(-l)=l,m(0)=2) = ( ^ 1 ] . (28) 



This operator obeys P^ 2 -j = /. That this is quite general, viz Pj — I, \fj is shown in [? ]. In appendix A we show 
that Pj = I \/j implies the operator relation, 

A a A' a = ^ A a . 

It is, perhaps, of interest that, if we associate the CB states with the position variable, q, of the continuous problem 
and its Fourier transform state, viz b=0 (cf. Eg. (1101) ). with the momentum, p, we have that the line of the finite 
dimension problem is parameterized with "initial" values of "q" and "p" i.e. m(-l) and m(0). 



VI. MAPPING ONTO PHASE SPACE 



We now define a mapping of Hilbert space operators, e.g. an arbitrary operator, B, onto the phase space - like lines 
of DAPG. The mapping is defined by [27j . 

B => V(j; B) = trBP-j. (29) 

Here Pj is a line operator within DAPG. (Alternatively we could have cast the mappings within APG as is clear from 
the discussion in the previous section.) The density operator may be expressed in terms of V(j; p): 



V(J; p) is quasi distribution [27j in the phase space like lines and points of DAPG. Thus, for example, the expectation 
value of an arbitrary operator B we have, 



1 d 

trpB=-J2v(j;p)V{j;B). (31) 



The quasi distribution may be reconstructed from the expectation values of the point operator A a i.e. MUB state 
projector's expectation value (obtained, e.g. by measurements), 



d 

Thence 



tr P A a = -J2 V(j; p)V(j; A a ) = -^2 V(j; p)K a>j . (32) 



\ E Y. V ^P)^ = V(j;p). (33) 
These equations are the finite dimensional Radon transform and its inverse. 

To summarize we give a comparative list of the central formulae: the expressions for the continuous phase space 
vs. the finite plane geometry ones: 



a. Phase space mapping, 



p^ W p (q,p), 
V p (j). 
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b. MUB state map (a=(m,b)) 



W\ x> 0) = S(x -Cq- Sp), 



A a =>■ trA a Pj = A 



c. Radon transform (a=(m,b)) 

Tl[W P }(x,9) = J^W p (q,p)5(x-Cq-Sp), (34) 
1 rf2 

nv P ](j) = jJ2 v p(j)^ tj . (35) 

Thus mapping of finite dimensional Hilbert space operators onto phase space-like lines and points of finite geometry 
was used to define a finite dimensional phase space Physics and applied to give the Radon transform and its inversion. 

VII. SUMMARY AND CONCLUDING REMARKS 

Finite geometry stipulates interrelations among lines and points. The stipulations for the (finite) dual affine plane 
geometry (DAPG) was shown to conveniently accommodate association of geometric lines and points with projectors 
of states of mutual unbiased bases (MUB). The latter act in a (finite dimensional, d) Hilbert space. This underpinning 
of Hilbert space operator with DAPG reveal some novel inter operators relations. Noteworthy among these are Hilbert 
space operators, Pj, j = 1,2, ...d 2 , which are underpinned with DAPG lines, Ly. it abides by P 2 — I Vj, and are 

mutually orthogonal, trPjPji = dSjji. These allow their utilization for general mapping of Hilbert space operators 
onto the phase space like lines and points of DAPG in close analogy with the mappings within the continuum of Hilbert 
space operators onto phase space via the well known Wigner function [HI. [33| . The physics of phase space involves 
Weyl-Wigner (W-W) mapping of Hilbert space operators to functions in phase space. The W-W image for p(x, 0) 
is the Radon transform of the Wigner function, W p (q,p). This comes about because the W-W map of the mutual 
unbiased state (MUB) involved, \x, 9), is S(x — Cq — Sp) which specifies the phase space point, q,p, that lie on the line 
stipulated within the S function. Thence the state reconstruction, gotten via the inversion of the Radon transform 
requires the properly normalized sum (integral) of all the lines going through that phase space point. This involves 
the handling of singularity in the continuum problem. Finite dimensional mapping onto the phase space-like lines of 
the dual affine plane geometry (DAPG) involves, likewise, mapping of an MUB state, \m,b). This map now relates 
DAPG points a = {m,b) and a lines ,j, via A QJ (Eq.(l2"3"|)') which replaces the continuum 5 function. Mapping and 
inversion in the finite dimensional analysis, in close analogy to the Wigner function of the continuum, is conveniently 
given in terms of a quasi distribution, V p (j) - the image of the state p. 

The approach allows an easy transcription from dual affine plane geometry (DAPG) to affine plane geometry (APG) 
and vice versa. 

Appendix: "Fluctuation Distillation" Formula 

Given, Eq(P|), P 3 = J2 a£j A a - I and, H3, i 3 / = /, implies 

a£j a'Ej 



Thus, 



Recalling that, Eq([T3]), A 2 a = A a allows 



E A a A a — 2 ^ A a . 



E A a A' a — ^ A a . 
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QED 
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